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ABSTRACT 


An analysis has been made of the corpressional wave spectrum 
excited in a fissioning uranium plasma confined in a cavity such as a 
gas cored nuclear reactor. Computer results are presented that solve 
the fluid equations for this problem including the effects of spatial 
gradients, nonlinearities, and neutron density gradients in the reactor. 
Typically the asymptotic fluctuation level for the plasma pressure is 
of order 1% . 
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1 . 


INTH3DDCTI0N 


The objective of this study has been to analyze the role of 

spatial gradients on the nonlinear development of fission power density 

driven instabilities in uranium plasmas. The principal impetus for 

study of such uranium plasma reactors has derived from the concept of a 

(1-5) 

high specific impulse nuclear gas-cored rocket engine. However, 

fissioning plasma reactors may also lead to development of a new class 
of magnetohydrodynamic generators ^ ^ , and perhaps may even have 
application to the direct nuclear excitation of lasers. 

The equations describing a uranium plasma involve the coupling 

of the plasma "fluid" dynamics to the neutron transport and also to the 

(12) 

radiation transport an the system . In such a reactor, the plasma 
would be enclosed (or nearly enclosed) by a neutron moderator — which 
is important in that it returns the thermalized neutrons back to the 
plasma core — while radiation transport dominates the transport of 
energy within the plasma. The system of equations used in this analysis 
is set out in the next section. 

The major application for the fissioning uranium plasma reactor 
(on which significant effort has been expended) is that of the gas- cored 
nuclear rocket engine shown belcw. In general, this system consists 
of a chamber in which a fissioning uranium plasma (at "5 ♦ 1Q 4 o K) is 
maintained in the central region and radiatively heats the propellant 
hydrogen gas which flews around the main core and is subsequently ejected 
from the rocket nozzle. The chamber is surrounded by a moderator which 
slaws dawn the fission neutrons and returns them to the -uranium core. 

/ 



50,000°K 



Figure 1.1 

It' is clear that if instabilities occur in this device, leading 
to either a fine-scale plasma turbulence or larger scale oscillations, 
they could well have undesirable results — such as enhancing the uranium 
fuel ejection from the rocket. An understanding of these processes (and 
their possible controllability) is of basic importance for the feasibility 
of constructing such plasma-cored rocket engines. 

In a mere general context, it should also be noted that the 
direct excitation of plasma motion and plasma waves by a fissioning plasma 
reactor may have useful applications to other energy conversion devices 
such as MHD generators. 

One source of instability is the growth of sound waves in 
the plasma due to the fission pewer density, P^ sg . The instability 
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medhanisra involved can be understood as follows: Consider a standing sound 

wave in a bounded region of fissioning plasma with a constant background 
thermal neutron density. In the wave compressions the fission power den- 
sity, P fj _ ss r increases due to the increased uranium density, and in the 

rarefactions P... decreases. This results in an increased pressure grad- 
fiss 

* T 

ient associated with the wave which in turn leads to a transfer of fission 
power to -the wave. It occurs because the wave compression tends to expand 
more rapidly than it was compressed. However, competing with this is the 
fact that radiation tends to transport the extra thermal energy out of the 
wave compressions . Radiation diffusion smooths out the temperature fluctu- 
ations of waves more rapidly the shorter their wavelength. This results in 
a critical wavelength, ^ cr ^ t / below which waves are stable and above which 
they are unstable. 


Last year we calculated the nonlinear evolution of such unstable 
waves to determine their limiting amplitude. The mechanism that finally 
limits their amplitude is the mode coupling of unstable waves in the wave- 
length range X > to radiation daitped waves in the range X < 

However in order to keep the equations tractable it was assumed that the 
uranium plasma was spatially homogeneous, and that only one mode was 
unstable. 


The work reported here relaxes these assumptions and deals with 
the full system of equations including the physical effects of nonuniform 
plasma density, temperature , neutron density, and radiation flux, on acoustic 
instabilities in the reactor. The teirperature dependence of the Rosseland 
diffusion coefficient and internal ionization energy storage in the plasma 
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was also included. 


The system of fluid equations describing a U-plasma heated 
by fission and losing energy by radiation diffusion is given by Eqs. 

(2.1) -(2.3) in the next section. Our procedure is to first choose a 
stationary solution to these equations corresponding to a given steady- 
state reactor. Then, a computer solution is obtained which follows the 
evolution of a perturbation in the plasma. 

Our general conclusion is that for typical parameters (of 
interest for gas core rocket engines) , the pressure fluctuation level 
due to fission-driven acoustic instabilities is of order 1%. The fact 
that this is considerably larger than that calculated last year can be 
assigned to the capability of the computer model to include the effects 
of spatial inhomogeneities for a reactor in which several modes are also 
unstable. These complex situations could not be treated by analytic 
methods . 

The above fluctuation level should not present a severe problem 
for open cycle reactors. Further, under some conditions (see section 5) , 
the partial transmission of sound waves out of the reactor cavity would 
have an important stabilizing effect which could suppress the fluctuations 
to belcw the 1% level. In an open cycle reactor, acoustic compression 
waves can propagate out of the cavity along the gas column ejected through 
the rocket nozzle (and a small amount of transmission can also occur 
through the walls) . This competes with the fission power density effect 
that tends' to energize such waves . 


"I 
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2 . FLUID EQUATIONS jjQR URANIUM PLASMA 


It is useful for conputational purposes to write the system 

of fluid equations for a U-plasma in conservation law form. We start 

... ( 12 ) 

with , 
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at 


+ V 




( 2 . 1 ) 


3V 
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3t 






MN 
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[N u KT u 


(1 + Z)] = 0 


( 2 . 2 ) 




I "u^u (1 + Z > + V 


+ J (i + Z) lycr^ • V u 


+ \£l = P f iss + I ’ ^ 


( 2 . 3 ) 


where the fission power density/ P fj_ ss ' is given by 


fiss 


N W v a Q 
u n o 


and 


N 


u 


V. 


u 


Z 

M 

T 


u 


N 


n 


v 


= nurtfoer density of U atoms and ions 
= radial fluid velocity 

= mean charge on U atoms or ions in units of e 
= mass of U-atom 
= U-plasma temperature 
= neutron number density 
= neutron thermal speed 
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Q = energy per fission event 

a = induced fission cross-section 

= radiation diffusion coef . (s a thermal conductivity) 

& = average energy used to ionize the U atoms to their 
average charge state Z 

The quantities £ and Z are shown ^ ^ as functions of T and pressure 
in figures 2.1 and 2.2. 

The last two equations can be re-written 

|t (MNV^ + V • [MJJ/Vu + £ (1 + Z)] = 0 (2.4) 


} =P £iss' <2 - 5) 

thereby expressing their local conservation-law form. 

Considering for example the motion of a plasma confined to 
a slab of thickness L with rigid walls, L.o ,. , 

V u (x = 0,L) = 0 , 

the instantaneous heat flux out of the reactor is 

ergs/sec/cm 2 . 

x=0,L 


Tot 





| x lJ-MNV 2 +|-NKr (1 + 
3t 2 u u 2 u u 


Z) 


+ N £ 
u 


+ 


V • I V J MSI V 2 + |- V N KT (1 + Z) 
— j -u 2 uu 2-uuu 


N V £ 
u-u 
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This is a radiation flux and is assumed to be carried off at the outer 
wall by for exairple a flew of cooler absorbing fluid. We will be inter- 
ested in an initial value problem in which the system starts at t = 0 
with small perturbation about a steady state solution. A spectrum of 
waves is then expected to be generated which reaches a steady-spectrum 
state. 


Equations in Dimensionless Form 

The following dimensionless variables are useful: 
x 



T = t/(L/V Q ) 


( 2 . 6 ) 


vhere L is a length in the system (e.g., sphere radius or slab thick- 
ness) , and 


V 


-KT o {1 + Z Q ) 


M 


P 

. o 


where P . p , T , Z . N , are all evaluated at the central point in 
o o o o o 

■the system. Thus we define 
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Further, define 


e = 


e = 


1 + z 
1 + z 


KT 0 (1 + Z c > 


( 2 . 8 ) 


In terms of these variables the system becomes 


9 

37 


P 

PV 


\ / 

3 


\ j PV 2 + j pT0 + pe ) 


9X 


PV 

pW + I pT0 


\ 


pV 2 V + | VpT9 + pVs - IJ, 

\2 2 . V 3 L dX 

o o 


0 
0 

L p fls S /Po v o 


(2.9) 


The dependence of the radiation diffusion coefficient and fission power 
density on T, N is given below. 


Radiation Diffusion Coefficient 


(12-16) 


The radiation diffusion coefficient is given by 


_ 16 a B T u 

^ - 3 k* 


( 2 . 10 ) 


where a is the black body constant and the Rosseland mean opacity 
in cm -1 . 
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Now there is sane uncertainty concerning the T u and pressure 
dependence of k^. Sane theoretical values are shown below, which accord- 
ing to Patch are expected to be correct to within an order of magnitude. 
Typically we shall be interested in the teirperature range from 50,000 
°K (reactor core) to 5000 °K (reactor edge) . In this range we see that 
k^ goes through a maximum. In view of uncertainties in the value of 
, and the fact that the T 3 dependence in (2.10) is stronger than the 
k R (T) dependence , we shall simply assume k^ to be constant in the follow- 
ing calculations. It would not be difficult to relax this assumption. 
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The quantity occurring in the dimensionless equation (2.9) 
involving radiation diffusion has a coefficient 


*R T o 

p o V o L 


16 a B T o 
3k R 


p v?l 

k O O 


T 


3 - 


KrT 3 


( 2 . 11 ) 


where we have defined a dimensionless diffusion coefficient K^. 
This can also be written 


T 


^ 


Kro 


LpV” 
o o 


V N K(1 + Z)L 
o o 


( 2 . 12 ) 


where is evaluated at the reference point x Q where T u 


V 


N u = N o • 


Fission Power Density 


The dimensionless quantity involving fission in (2.9) is 


lP V 3 
o 


\ T> — I 

' L 

1 riss \ 

^ Po V o ) 


N N v crQ 
u n o 


For the neutron density in the cavity we allow a space dependence and 
write 


N =7) g (X) 
n n ^ 


where g (X Q ) = 1 evaluated at the reference point x Q . Thus, 
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p V 3 
H o o 


P-. 

fiss 


where the dimensionless fission power is 


P 

riss 


P V° 
K o o 


v o°Q R? 7n 

3 flss '° KT (1 + Z )V 
o o o o 


P f iss ,o 


Wo 00 


is the fission power at the central point x Q where = N q and 
T u = T o‘ 


Equations for Slab or Sphere 


In surarnary then, equation (2.9) can be written as 


j pV 2 + j pT0 + pe 


pv V + I pT0 
J pv 2 v + | VpT0 + pVe 


rr m3 8T 

3X 


P fiss P^ (X) 


This equation is simplest to integrate for one-dimensional geometries 
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such as slabs, spheres, or cylinders. It becomes complicated for two 
dimensions, and only a few 3-D fluid codes have been written to date and 
invariably consume large quantities of machine time or the use of extremely 
coarse grids. Two one-dimensional cases of interest to us are: 

Slab 

/ pv \ 

pv 2 + pTe 

^ P V S + §V P T8 + PVe-K^H j 

(2.16a) 

/ 0 \ 

0 

\ fiss 13 ^ / 


3 

3x 


\ 


\ 


PV 


9 

9X 


jpv ! 


+ pT9 + pe 


/ 


Sphere 


3 

3t 


/ 


R 2 p 


pv 

J pv 2 + | pT0 + pe. 


3x 




R 


/ 


0 \ 

0 


' R2 ^fiss p 9 1 


R 2 pV \ 

pV 2 + pT9 

I P v 3 + | vpTQ + pVe - i^T 3 || / 

(2.16b) 
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It should be noted that the quantities Z and £ are also 
functions of T and P. They occur in (2.16) in the dimensionless 
variables 


8 - 1 + 2 £ _ £ 

1 + Z Q ' & KT o (l + Z Q ) 


(2.17) 


Now the temperature dependence of Z, £ is more rapid than their pressure 
dependence. Also, the reactor plasma will have a pressure which is 
nearly constant as a function of position (except for fluctuations 
associated with wave motion) , whereas the plasma temperature varies 
by an order of magnitude from the center of the reactor edge. Thus it 
is reasonable to set 


Z(T, P) £ Z(T, P ) 
£(T, P) = £( T, P Q ) 


(2.18) 


in the fluid equations, i.e., only account for the temperature depend- 
ence of Z, £ in the time-dependent problem. 

This has been done by making a linear fit of the UPLAZ-2 
(4) 

computer calculations for these functions, i.e., 

Z = a + b T = a + b T T 
u o 

(2.19) 

S — c+dT^ = c + dT Q T 

The coefficients a, b, c, d are functions of pressure evaluated at the 
approximately constant reactor pressure. 
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3. STATIONARY SOLUTIONS 


From equation (2.16) we see that the stationary solutions in- 
volve solving 


pT0 = 1 


(3.1) 


and 


w i- 

*R 3X 


T 3 

1 3X / 


+ P fiss p ^ (X) 


= 0 (slab) 


(3.2) 


or 


*R 


W ( t!r 2 H) + p fiss p 9(WR 


2 _ 


= 0 (sphere) 


(3.3) 


Using (3.1), these become 


(1 + Z)T 


3X 


K) 


+ 5 g(x) = o 


where 


5 = 


P-. (1 + Z ) 

f xss o 


i.e. , 


*R 


(slab) 


(3.4) 


(1 + Z)T ^T 3 R 2 || j + CR 2 g(R) = 0 (sphere) (3.5) 


/ 3T \ 2 

W I 


+ t — + -£-a = 0 


3X 2 (1 + Z)T 3 


(slab) 


(3.6) 


(ir-s 


+ + ~ = 0 (sphere) 

9R , R 3R (1 + Z ) T 3 


(3.7) 
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Now writing the grid spacing as h, it follows that 


3T 
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3 2 T 
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+ 



(3.8) 


Thus (3.6) or (3.7) give 


a quadratic in T j+1 


I 

where 


slab 
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T j + 1 - k ( - B + (B 2 - 4SC) * ) 


(3.9) 
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2 T 2 -2T 2 + T.T. i K h 2 g (X = jh) 
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4 
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j 3 


3 t 2 
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3 3 3-1 


T m , K h 2 g(R 4 ) 
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3 o 3' 
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where we have noted equation (2.19) for Z . 

Sane sample curves are shown in Figs. 3.1 to 3.4 for a slab of 
thickness L (L is taken as the reference length in Eq. 2.6). The refer- 
ence point where N = N q , T^ = T q is taken as the center, i.e. , x = L/2, 
X = . 5 , and the neutron density function g was taken to be 

g = e “ <X 'I )2 5e CL 2 (X- .5) 2 _ 


The constants used in Figs. 3.1 - 3.4 are: 


a = - .4375 
b = 9.37 • 10" 5 


Z = 4.247 
o 

T = 5- 10 4 °K 
o 

OiL 2 = 2 or 0 


N = 1.5 10 19 
o 

k^ = 1.2 10 2 on -1 


c = - 3.17 • 10 
e = 2.64 • 10“ 15 


-1 1 


} 


for Z = a + bT T 
o 


(Scale for neutron density) 


for £ = c + dT T 
o 


Various values of £ which determines P^. „ „ . 

f iss ,o 


= 3.2 • 10 8 er/sec/deg. 


(3.12) 
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4. TIME-DEPENDENT PROBLEM 

The program that integrates (2.16a) makes use of the Lax Wendrof 
two step method. Equation (2.16a) has the form 


It “ 2 


(4.1) 


Step 1 : 

Compute provisional values at centers of rectangular meshes of 
the net in the x - x plane using 


- k <A, + (ft - & (F 1 ? . - f?) + 21 (rf? + o'?) 

~3+h 2 — g+l —j 2Ax — d+1 —j 4 -g+1 —j 


(4.2) 


Step 2 ; 


Then advance in time using 


u ^ 1 

-3 


s? - i ®j+>s + 2j->s * • 


(4.3) 


The walls are treated as rigid so that the fluid flux pV 
vanishes there. However, the energy flux is a floating quantity at the 
wall, i.e. , it is continuous across the boundary. A spatial Fourier 
analysis is also made of the various fluid variables, 



E (A sin k x + B cos k x) 
n n n n 

Trn 

L 


(4.4) 
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Interpretation of Fourier Mode Behavior and 
Comments on Turbulence Equations 


When the stationary system described in Section 3 is perturbed, 

a spectrum of sound waves is generated. Ibis spectrum settles down into 

an asymptotic finite amplitude spectrum after about 50 sound transit times 

across the reactor. It derives its energy frcrn the fission power driven 

(12) 

instability described earlier. 

Now the presence of such a wave spectrum also causes the time- 
average tarperature and density profiles across the reactor (about which 
fluctuations occur) to depart from the steady state solutions of section 3. 
One can clarify this as follows: 


First define a time average over a time scale larger than the 
wave periods. 


<> 


1 

T 


T/2 


f dt 

-T/2 


(4.5) 


with T » w _1 . 

n 


The averaged fluid equations (2.16) then become 


9T 
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P 
pV 

i pV 2 + | pTG + pe 


> + 


9 

9X 


< 


PV 

pw + IpT0 
2 pv 2 v + | vp T0 + 


+ pve - Kj, T S § I > 


= < _ 0 > 


P £iss pg(X) 


(4.6) 
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This equation, together with (2.16), can be used to generate equations 
for the averages <p>, <V>, <T>, etc., and fluctuations Sp, 6V, ST, etc. 

In the 'fetationary turbulence" state (and for one-dimension) it 

follows, 

< Sp <5V > = 0 

<p><T><0> + <p>< 5V 2 > + < Sp ST > < 0 > + <(5Ta0><p> 
+ < Sp 50 > < T > + < Sp ST 50 > = <P> = constant. 


■[ j < Sp 6V 3 > + | < p > < SV 3 > + | < Vp T0 > 

+ < pve > - K* < T 3 f >} = < P f . ss pg > 


(4.7) 


We now see that the stationary solution to the above system differs fran 
the stationary solution to (2.16) because of the various pressures, fluxes, 
etc., deriving from the fluctuations. 

Returning to the exact system solved by computer, the fluid 
variables are expanded 


P?V,T = L 


{^ P sn' 


V , 

sn' 


T ) 
sn 


sin 


k x 
n 


(p cn' 


V , 

err 


T cA 


cos 


k x 
n 


} 


(4.8) 


where p , V , T are functions of the "slow" time variable T and 
w oc oc oc 

all other p , V , T are functions of both slow T and also oscillate on 
n n n 

the fast t-scale. Also, 


< p, V, T > = ( p Q , 0, T q ) 


(4.9) 
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Pressure Fluctuations 


Consider the fluctuating part of the pressure (including wave 


pressure) , 


<5P = P - < P > 

= pV 2 + pT0 - < pT 0 + pV 2 > 

co (4.10) 

5 E 5P n 

n = 1 


The Fourier canponents 6P n all fluctuate about zero, as illustrated below: 



They are particularly simple to evaluate at the wall where V = 0 and 


oo - 

®w - T w 0 w L { p ns sln k n x + p ne 003 k n x } I 4 ' 11 ) 

n=l -* 
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An alternate quantity of interest (which we will in fact plot 
using the oornputer results) , is the total pressure fluctuation about the 
initial state , i.e., 

6P = pV 2 + pT0 - (pT0+ PV 2 )^ (4.12) 

In this case the various spatial Fourier components of 6P oscillate in 
time about an average pressure that drifts away fran the initial pressure 
due to the wave pressure terms, i.e., Fourier components appear as 



Figure 4.2 

These are the quantities that will be plotted for various neutron gradients 
and reactor parameters in Figs. 4.2 - 4.7. They are graphed for the one- 
dimensional slab case. The asymptotic Fourier amplitudes are all about 
1% for the cases considered with an enhancement of approximately a factor 
of 2 for the cases in -which the neutron density gradient was nonuniform. 
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T = t/(L/V 0 ) 

Figure 4.5 

Evolution of modes n = 1, 2, 3 for E, = 1.2, ctL' 
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Logarithmic, 2x2 Cycles 












T = t/(L/V Q ) 

Figure 4.6 

Evolution of modes n = 1,2,3 for £ = 1.4, aL' 
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Logarithmic, 2x2 Cyck-s 




































Relation, between Spherical and Slab Geometry for the 
Fissioning Plasma Problem 


lated. 

Slab 


The equations for spherical or slab plasmas can be simply re- 


JL 

3x 


P 

PV 

| pV 2 + | pT0 + pe 


ax 


/ 


\ 


pv 

pV 2 + pT0 
I pv 3 + | VpT0 + pVe - KpT 3 || 


\ 


0 
0 

P fiss ^ 


\ 


(4.13) 


Sphere 



( r2 p \ 

a 

( R 2 pV ' 

a 

3x 

pv 

i i 

+ 8R 

pV 2 + pT0 


\ R 2 [ | pv 2 + | pT0 + pe ] ' 


^ R 2 [ j pV 3 + | VpT0 + pVe- K^T || ] / 



/ „ \ 


(4.14) 


\ R2p fiss ' 


V7here p, V, T, e, K^, P fj_ ss ' axe the dimensionless density, fluid velocity, 
tenperature, internal ionization energy, radiation diffusion coefficient. 
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and fission power density. The function g(X) is the neutron density 
in the cavity. 

Now in (4.14) consider the variable change from R, V to X, W 

where 


X = R 3 , dX = 4 ttR 2 c3R 
W = 4ttR 2 V 
In terms of X, W, equation (4.14) be cones 


(4.15) 



(4.16) 


This equation is now identical with (4.13) [where W now takes the place 
of V ] , except for the four underlined terms. Three of them replace 
occurrences of V 2 in the energy density with calculations from W; the 
fourth is a temperature difference in R instead of X. Thus, to solve the 
spherical problem with the slab solution cne must make the replacements 



- 34 - 




where underlined 


. Note, 


_W = W 

4ttR 2 (4rr) V 3 (3x) 2 / 3 


(4.17) 


These relationships are useful in developing hydro-codes for the slab 
and sphere problems. 


Garments an Radiation Flux Instabilities 

It has been suggested that under sane circumstances radiation 

(17) 

flux may drive acoustic instabilities in a U-plasma (see for example 
"An Acoustic Instability Driven by Absorption of Radiation in Gases" by 
M. J. Monsler, MIT Tech. Note CSKT-69-4 (1969)') . Our set of equations 
(2.16) contains this physical possibility since a .large radiation flux is 
of course present in the steady state reactor solution that is perturbed 
at t = 0. 


It is fairly difficult to make a reliable analytic decision 
concerning this possibility, and Monsler 1 s analysis had several weaknesses. 
However, in order to extract seme limited information on this we considered 
a simple model for sound wave propagation in a medium that is treated as 
homogeneous except in terms involving the radiation flux. The linear 
equations are then 



N 


o 



(4.18) 







0 


(4.19) 


35 



to|w 


(4.20) 


K 


f 3T i 9N l ^ 

( N o TE + T o at" y 


, 3V, a / 3 (T + T ) \ 

+ fw 3ST = ^K-fx — J 


- a 


3T X 

~3x 


+ 3 T i + Kr 


3 2 T. 


3x 2 


It follows that modes e 1 ^ ^ have a dispersion relation 


a)(3w 2 - 5kV 2 ) = - (w 2 - kV 2 ) (-ika + 3 - R^k 2 ) 


(4.21) 


where v 2 = KT^/M. Assuming radiation has a small effect one can write 

to = w + w. and treat w. as small with w = /5/3 kv. This gives 
O' 1 1 o 

(3 - Kr)c ) 

o 


The range of stable wavenumbers is 


k 2 


crit 


1_ 

V 


3 2 K 


R 


3T 2 

o 



(4.22) 


X«0* f 

v ~ /n(n - 1) 

K crit \ 

where & T is the gradient scale of T^ and it has been assumed K^a T 11 . 

Now the approximations used in the analysis break down except 
for the wave number range k > several k cr j_ t / i»e. (similar to a WKB 
assumption) , short wavelength acoustic waves are stable . The analytic 


- 36 - 



question of the radiation stability of long wavelength waves with wave- 
lengths however remains open. Our computer studies bear out the 
conclusion that the radiation flux does not appear to give rise to in- 
stability. 



5. STABILIZING EFFECT OF SOUND WAVE TRANSMISSION 
OUT OF THE REACTING PLASMA 

Instability occurs because the average effect of over- 

comes the damping effect of radiation diffusion for sound waves in the 
wavelength range \ > However, the growth rates are small (typically 

10“ 2 w) , so that it is of interest to consider the damping effects due to 
escape of sound energy fron the fissioning plasma. 

In realistic situations sound wave energy is partially trans- 
mitted through the container walls (as well as down the orifice gas column 
in the rocket case) . The question arises of how to represent this trans- 
mission coefficient in the equations. 

New the equations involve the fluxes and energy flux F. 

Obviously the wall boundary condition N^V^ = 0 must not be changed since 
mass is conserved in the reactor even though energy is partially transmitted. 
The boundary condition representing transmission can most easily be intro- 
duced into the energy equation as a sink of sound energy density at the 
wall, i.e., rewrite the energy equation as 




Various models for S are possible. For example, a simple non- 
dispersive one is 


S = a y t) - E o 


(5.4) 


where E q is the energy density for the stationary solution (taken as 
the zero-frequency part of E ) , and (t) the finite-frequency com- 
ponents of E evaluated at the wall. In this case if one considers the 
simple limit P f£ ss = 0 and = 0 , (5.1) becomes 




6(x - x ) 
w 


Since (as has already been assumed) , E is nearly spatially homogeneous, 

r L . 

an integration / dx removes the 6-function and gives 

J o * 


If = - E < E - e o ) • «•« 

Note the wall boundary conditions on which the code is based lets heat 
out of the system but not back in. Thus we consider only perturbations 
for which E > E q . Equation (5.5) becomes 

_ a t 

E - E = (E (t = 0) - E )e L (5.6) 

o o 

i.e., the system decays to its original stationary energy state E q due 
to the escape of sound wave energy in a time L/a . 

For the more general problem we will use (5.4) in equation (5.1). 
However it is useful to relate a to the transmission coefficient. 
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Figure 5.1 


Assume a fraction A of the energy of a traveling sound wave 
is transmitted as shewn, and that A is independent of frequency. In 
the perfectly reflecting case A = 0 , the standing wave solutions for 
the slab have amplitudes that vary in time only due to and radi- 

ation diffusion. For A > 0 however, their amplitude decreases due to 
escape of their decomposed traveling waves. The strength of the sink S 
in (5.1) can be written 


S = A | (E - E ) V (5.7) 

2 W OS 

where V is the sound speed at the wall and (E - E )/2 is the energy 

S W O 

density in the right-traveling wave. Thus the exponential decay rate in 
(5.6) becomes 

)tl ■ . 
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Critical Transmission of Coefficient for Stabili 


Transmission of sound energy out of the reactor cavity will over- 
coma its amplitude growth inside the cavity if 


< or 


(5.8) 


where is the linear growth rate. This relation leads to a critical 
transmission coefficient such that if A > A cr _j_ t the system is stable, i.e. , 


A > A 

cret 


A < A . , 
cret 


(stable) 

(unstable) 


(5.9) 


It should be pointed out that this relation is fairly accurate, but the 
exact computer solution with a generalized boundary condition allowing 
transmission will yield a more precise value for A^^ .’ 

Now for reactors of interest, - kV g , so that the 
appropriate growth rate calculated earlier is 


i(k crit ‘ k )K Ro 


(1 + Z o + N oV 


(5.10) 


3N K 1 + Z + T. Z + — 
o L o o T K 


Also, several modes are unstable with the fastest growing mode correspond- 
ing to k a minimum, i.e., we can write k 2 .. - k 2 - k 2 and note 

cret cret 
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R U + z o + T o y 


k .. = 
crit 


P fiss/Q 


r KT 

K V 2 1 - (1 + Z + N ZJ 

itos V 2 M o o N 

R 


(Notation is as in our earlier report, ' ' or the 1971 final report) . 


71 P fiss,o 
w = — 

3MN V 2 
o s 


n = 


1 +Z o +T c/t 

1 + Z +TZ m +^- 
o o T K 


It follows that 


crit 


. = (^) 
it \V J 


P,. 

fiss,o 

3MN V 2 
o s 


or, since = kV s = ttV s /L (for the slab) , 


A = Tr- 
ent l UL. 


1 / fiss,o 

\ 1 3MN V 2 
V os 
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Typical values are 


p £ . = 1.3 * 10 10 ergs/am 3 /sec 

fiss,o 

MN = 6 • 10 -3 gm/ati 3 

o 

V = 2.7 • 10 5 an/sec 
s 

L - 10 2 cm 


so that ~ ^ • 10 3 . 

For example, consider the open cycle reactor in which waves are 
trapped in the cavity except when traveling out the nozzle as shewn. An 
approximate effective transmission coefficient is 



Figure 5.2 

i.e., we would expect the escape of sound wave energy to stabilize the 
acoustic mode provided A^^. < ^ e ff ec tive ' ^"* e ‘ ' r nozzle > R cavity. 
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